The geometrical nature of optical resonances in nanoparticles 
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We give a geometrical tiieory of resonances in Maxwell's equations that generalizes Mie formulae 
for spheres to any dielectric or metallic particle without sharp edges. We show that the electromag- 
netic response of a particle is given by a set of modes of internal and scattered fields and reveal a 
strong analogy between resonances in nanoparticles and excess noise in unstable macroscopic cavi- 
ties. We give examples of two types of optical resonances: those in which a single pair of internal 
and scattered modes become strongly aligned in the sense defined in this paper, and those resulting 
from constructive interference of many pairs of weakly aligned modes, an effect relevant for sensing. 
We demonstrate that modes can be either bright or dark depending on the incident field and give 
examples of how the excitation can be optimized. Finally, we apply this theory to gold particles 
with shapes often used in experiments. 
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The interaction of light with wavelength sized particles 
has been intensely investigated for more than a century, 
continuing to provide interesting and surprising results. 
This coupling is essential in single molecule spectroscopy 
and photon processes [H [2] , while interference between 
different scattering channels of particles produce classi- 
cal analogues of quantum processes [31 3] , and carefully 
designed particles underpin the physical realization of 
metamaterials [5H3- The basis all of these effects is the 
strong variation of particle-light interaction with wave- 
length. For particles much larger than the wavelength of 
light, resonances are described by closed orbits of light 
rays |5] inside the particle. This geometric approach be- 
comes less and less effective as the size of the particle 
decreases, eventually requiring the solution of Maxwell's 
equations. Mie- type theories [51 [TO], based on symmetry 
and coordinate separability, provide analytical descrip- 
tion of resonances for a few specific shapes of particle. For 
spheres, a resonance occurs when the coefficient of one 
of the electric or magnetic multipoles in the field expan- 
sion becomes infinite for particular values of the particle 
radius, permittivity and susceptibility. In this frame- 
work, resonances are not connected to anything similar 
to energy levels (eigenvalues of the familiar hermitian 
operators of quantum mechanics), and do not show an 
obvious geometrical interpretation. Here we introduce 
an hermetian operator to reveal the geometrical nature 
of resonances in the Maxwell equations. 

We consider metallic and dielectric particles without 
sharp edges and of dimensions of the order of the wave- 
length of the incident field or smaller, so that the in- 
teraction between light and matter inside the particle 
is described by local macroscopic permittivity and sus- 
ceptibility. The tangential components of electric and 
magnetic fields E, H are continuous on passing through 
particle boundaries, and the energy scattered by a par- 
ticle flows towards infinity. We use [H] six component 
vectors F = [E, H]'^ for electromagnetic fields; their pro- 



jections, /, onto the surface of the particle are surface 
fields that form a space T-L where scalar products are de- 
fined in terms of surface integrals. In this formalism the 
boundary conditions become 

f = f-r, (1) 

which has a simple geometrical meaning in H: the pro- 
jection /° of the incident field onto the surface is equal 
to the difference between the projections of the internal 
and scattered fields, /* and f. This suggests that an 
external field with small tangent components can excite 
large internal and scattered surface fields provided that 
these cancel each other. This happens when the "angle" 
between these two fields, and therefore their difference, 
is small. 

Importantly, for the particles considered here, sets of 
electric and magnetic multipoles that are linearly inde- 
pendent and complete in H can always be found [121 IE] • 
Using these sets we then form separate orthonormal bases 
{in} and {sn} for the projections of internal and scat- 
tered fields. These fields, which we call principal internal 
and scattering modes, are orthogonal to all but at most 
one function in the other space. For spherical particles, 
each pair of modes corresponds to a pair of electric or 
magnetic multipoles of Mie theory. For non-spherical 
particles, principal modes are instead combinations of 
different multipoles (although in some cases there can be 
dominant contributions from a specific multipole). We 
can find the internal and scattered fields for a particu- 
lar incident field by decomposing this field into a sum 
of pairs of modes /° — X^n '^n*" ~ o,f^Sn, where a'/'' are 
amplitudes of the internal and scattered principal modes 
for that specific incident field. The angle, between s„ 
and in is defined as 

i„ • s„ = cos (^„), (2) 

where we choose the arbitrary phase factors so that the 
scalar product is either positive or null. The terms on the 
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right-hand side of eq. ^ are the principal cosines intro- 
duced for the first time by Jordan in 1875 [M]. sin (^„) is 
the orthogonal distance between s„ and and cos(^„) 
their statistical correlation [TS]. These angles are invari- 
ant under rotation and/or multiplication by a phase fac- 
tor [IS] and characterize the geometry of the subspaces 
of the internal and scattered solutions in H. This ge- 
ometry is induced by the particular particle through the 
surface integrals of the scalar products; its relevance to 
scattering and resonances has not been previously real- 
ized. The importance of the principal cosines is twofold: 
Theoretically they provide analytic equations for the co- 
efficients of the internal and scattered principal modes, 
generalizing the Mie formulae and clarifying the nature 
of all scattering channels of a particle. Numerically, they 
allow us to reduce large matrices to their sub-blocks and 
eliminate the need of numerical inversion for the deter- 
mination of the mode coefficients. 

Minimization of |/° — J^n'^n^n + aJi^nl leads to the 
equation 
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where X, S are matrices whose columns are formed by the 
principal modes {«„} and {s„}, respectively, 1 is the iden- 
tity and a*, a" are the coefficients of the principal modes 
in the field's expansion. The most important part of our 
theory is that principal modes are coupled pairwise, i.e. 
that the matrix C is diagonal. The interaction of parti- 
cles with light can then be interpreted in terms of eigen- 
values and orthogonal eigenvectors wi^ = («„ ± Sn)/\/^ 
of the hermitian operator in Eq.([3|. However, away 
from the surface one measures either internal or scat- 
tered fields, so we transform the eigenf unctions, {w^}, 
to find the coefficients of the principal modes: 
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Here 4 = z„ - cos(^„)s„, = s„ - cos(^„)i„ are 
bi-orthogonal to i„,s„ (i^ • s„ = s'^ ■ i„ = 0) with 
i'n ■ in = s'n ' = sin^ (^„). aj, , fl^ in Eqs. are found 

by projecting the incident field onto non-orthogonal 
vectors i„,Sn and sin(^„) is defined as the Peterman 
factor [17] that gives the order of magnitude of tran- 
sient gain and excess noise in unstable cavity modes. 
Therefore the presence of strongly aligned vectors with 
sin(^„) << 1 is at the origin of large surface fields in 
nanoparticles as well as large transient gain and excess 
noise [18] in macroscopic unstable cavities and dissipa- 
tive systems governed by non-hermitian operators |19j . 
However, the physical origin of the non-orthogonality of 
i„,s„ is not dissipation, but the fact that internal and 



scattered fields are not modes of the hermitian opera- 
tor and are correlated, even in the case of non-absorbing 
particles. 

These surface modes and their amplitudes fully specify 
the electromagnetic behavior of the particle everywhere 
in space, for example electric currents within the particle 
and its surface plasmons can be found using the fields 
and the Ohm equation. Eqs. (|4][5| themselves have sev- 
eral important consequences. The left terms in the scalar 
products depend exclusively on the particle, this allows 
us to disentangle its properties from those due to the spe- 
cific incident field. As an example, we now consider reso- 
nances. The principal modes and scalar products depend 
on the frequency dependent permittivity and susceptibil- 
ity of the particle, and therefore the principal angles ^„ 
change with the frequency of incident light. Internal and 
scattered coefficients diverge when the denominators of 
eqs. Qjs]) vanish. This happens when a pair of normal- 
ized internal and scattering modes are parallel. For a 
sphere the angular dependence of internal and scattered 
modes can be factored out and the condition i„ — Sn can 
be recast in terms of the amplitude of the electric and 
magnetic components giving the Mie resonance condi- 
tion. For physical particles, the linear independence and 
completeness of the principal modes makes perfect align- 
ment impossible. So, as with spherical particles [2U], ac- 
tual resonances correspond to minimum angles 7^ 0) 
of pairs in T-L. One important difference from spheri- 
cal particles exhibited by non-spherical structures is that 
the total flux of energy scattered or absorbed (integrals 
of the Pointing vectors over all directions) is given by 
the sum of principal mode contributions plus interference 
terms between modes, which are absent for spherical par- 
ticles. Hence for non-spherical particles, strong peaks in 
the efficiencies can be caused by constructive interference 
within a group of modes. Furthermore, Eqs. (|4j[5]) show 
that modes can be "dark" for specific incident fields. The 
largest internal surface field is produced by an incident 
surface field parallel to i'^: For well aligned modes, the 
coefficient |a*|, ja"! are of the same order while, for weak 
resonances, ja"! is much smaller than |a*|, leading qual- 
itatively different absorption and differential scattering 
cross sections (DSCS). Similar effects happen for an in- 
cident surface field parallel to s^, exchanging |a*| and 

Regarding numerical calculations, we remark that 
Y^n=i'^n^ri,Y^n=i^n^n couvcrgc to the cxact ficlds for 
iV — >■ 00 at any point inside and outside the parti- 
cle [IHEl], and \f'^ + Y^n=i O'n^n — CLnin] allows US to ver- 
ify the accuracy of the calculations for finite N JIB] . Most 
importantly, the form of Eqs.(4|5) remains unchanged as 
N — 00, even if s„ change [16] . 

We now consider light interaction with some specific 
gold particles, using a fitted dielectric function [55]. 
Fig. [1] shows calculated near and far field properties of a 
rounded nanodisc (height 20 nm and diameter 120 nm). 
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FIG. 1: (Color online) Calculated modes of the gold nanodisc. 
a) Scattering, absorption and extinction efficiencies, showing 
the strong resonance at 613 nm. b) The DSCS as a function 
of wavelength, illustrating the dipole nature of this particu- 
lar resonance, c) and d) Excitation paths of scattered and 
internal fields over the principal angle landscape, e) and f) 
Near field plots of the dipole resonance disc viewed parallel 
and perpendicular to the electric field of the axial light. 



FIG. 2: (Color online) Calculated optical modes of the 
rounded nanorod when illuminated axially. a) Scattering, ab- 
sorption and extinction efficiencies showing the 205 nm mode 
and weaker absorption peak at 486 nm. b) DSCS showing the 
strong forward scattering of light by the particle, c) and d) 
Scattered and internal excitation paths, showing the compos- 
ite nature of the modes. 
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illuminated axially by plane-wave light. We find a strong 
dipole resonance (fig.[I}i,b) at 613 nm, in agreement with 
a previous experimental and numerical investigation |23j . 
Figs. ^:d) show the landscape of principal mode angles, 
which are intrinsic properties of the particle, as a func- 
tion of mode index and wavelength of incident light. The 
"height" of this landscape is sin(^)~-'^, the largest value of 
eq. [ijand eq. [sjfor |/°| = 1, while the amplitudes of the 
internal and scattered principal modes due to the specific 
axial incident field are shown overlaid on top. These ex- 
citation paths demonstrate that the observed resonance 
is a Mie-like single principal cosine pair, which comes 
into maximum alignment (its smallest value of ^) at the 
resonance. Other strongly aligned dark modes are also 
present, but not excited by this particular incident field. 
The internal field also contains a second more weakly 
aligned excited mode, where its counterpart in the scat- 
tering field is not excited. This mode and of the reso- 
nant mode have quasi-degenerate principal angles on the 
shorter wavelength side of its peak; this strengthens the 
internal field and explains the appreciable asymmetry of 
the absorption efficiency. The near field around the parti- 
cle at resonance (fig. [l^,f ) shows that the resonance is an 



electric dipole. For non axial incidence, the main peak is 
smaller due to a weaker coupling with the resonant mode. 

Fig. [2^,b) shows the calculated optical efficiencies and 
DSCS of a 480 nm long rod with diameter 40 nm, illumi- 
nated axially by plane- wave light. There is a strong reso- 
nance at 205 nm and a weaker absorption peak at 486 nm 
(an extremely weak corresponding scattering feature is 
present at ^^515 nm). The DSCS demonstrates that this 
particle strongly scatters this incident light forward, par- 
ticularly at short wavelengths. The peak at 205 nm is 
not a single mode, but instead the excitation amplitudes 
of the principal modes (fig [2j:,d) show it to be due to 
constructive interference of contributions from a group 
of several weakly aligned principal mode pairs. Collec- 
tively these modes very efficiently extract energy from 
the incident field. Such multimode resonances are poten- 
tially useful for sensing applications because interference 
between different scattering channels leads to enhanced 
sensitivity to perturbations near to the scattering sur- 
faces. The absorption feature at 486 nm is also due to 
a group of mode pairs, in which only the internal modes 
are excited, all becoming more strongly aligned at around 
this resonance. Fig. |3] shows the same particle illumi- 
nated equatorially with an incident light polarization of 
45° with respect to its long axis. Analyzing other parti- 
cles with the same diameter but varying length, we find 
that the broad feature at around 200-450 nm, contain- 
ing structures similar to the composite modes of Fig. [2] 
is insensitive to the particle length. It also shows no 
clear hot or cold spots (see figjsjj). The sharp resonance 
at 676 nm shifts with the length of the rod and its sur- 
face field has the strong nodal local structure (Fig. [3t) 
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FIG. 3: (Color online) The particle from fig. [2] but illumi- 
nated from the size with an incident light polarization of 45° 
with respect to its long axis, a) Far field scattering efficiencies 
showing the presence of both a broad feature similar to Fig. [2] 
at around 200-450 nm and also a sharp resonance at 676 nm. 
b) The DSCS for equatorial illumination c) and d) Scattered 
and internal excitation paths; the principal angle landscape 
includes additional modes that cannot be excited by symme- 
try for axial incidence (cf. Fig. [2|. The sharp mode is again 
a single well aligned principal pair, the best aligned of three 
(and the only one excited by this field), e) and f) Near field 
for the 480 nm rod, shown at the broad feature (207 nm) and 
at the "waveguide" Mie-like mode at 676 nm. 



property of the particle, and the way an incident field 
couples to these modes. In particular, there can be res- 
onant modes that are essentially dark with respect to 
a given incident field. The amplitudes of the principal 
modes specify the field both inside and outside the parti- 
cle. This geometric approach is applicable to a large class 
of processes, where the interaction between the system 
and the environment is described by a set of functions 
complete at the boundary of the system, including other 
scattering processes such as of acoustic or electron waves, 
and coupling to optical cavities. 




~200 'too 

Wavelength (nm) 







3D, 




















Mode ind9x 




+4-30 1 


T 


+240 




200 400 

Wavolgngth fnm) 



0.1 



-4BtJ -240 D +240 +4B0 
Posilion {nm} 



-480 -240 +240 +4aa 
Position [nm] 



of a "waveguide" mode on the long axis. The excitation 
paths (Fig. [3j:,d) are shown including additional modes 
that cannot be excited by symmetry for axial incidence 
(cf. Fig. [2]). As with the disc, the sharp resonance is again 
a single well aligned principal pair, the best aligned (and 
the only one excited by this field) of three such waveguide 
modes that shift with the length of the particle. 

In summary, we provide a geometrical explanation for 
resonances that generalizes Mie formulas for spheres to 
all isotropic particles without sharp edges and reveal 
strong analogies between resonances in nanoparticles and 
excess noise in macroscopic cavities. We show that there 
are sharp resonances caused by strong alignment between 
one internal and one scattered modes and broad res- 
onances due to several pairs of internal and scattered 
modes with weak, but similar, alignment. Our equa- 
tions indicate that any observation depends both on the 
structure of the principal modes, which are an intrinsic 
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